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Abstract. We report on observations of the Crab pulsar with the MAGIC telescopes. Our data were 
taken in both monoscopic (> 25 GeV) and stereoscopic (> 50 GeV) observation modes. Two peaks were 
detected with both modes and phase-resolved energy spectra were calculated. By comparing with Fermi- 
LAT measurements, we find that the energy spectrum of the Crab pulsar does not follow a power law 
with an exponential cutoff, but has an additional hard component, extending up to at least 400 GeV. This 
suggests that the emission above 25 GeV is not dominated by curvature radiation, as suggested in the 
standard scenarios of the OG and SG models. 



1. Introduction 

The Crab pulsar is the compact object left over after a historic supernova explosion that occurred in the 
yeai" 1054 A.D. It is among the brightest known sources at GeV energies. However, despite numerous 
efforts, a spectral steepening made its detection above 10 GeV elusive until 2008. Consequently, the 
presently prefered Outer Gap (OG) H and Slot Gap (SG) Jl models predict the emission to be produced 
through curvature radiation, which implies an exponential cutoff at a few GeV. The new satellite-borne 
gamma-ray detector Fermi-LAT, which can measure the spectra of gamma-ray pulsars up to a few tens 
of GeV, became operational in August 2008. The spectra measured by Fermi-LAT can be described by a 
power law with an exponential cutoff, which supported the OG and the SG model. 

In 2008, the MAGIC-I telescope finally detected the Crab pulsar above 25 GeV with a newly 
implemented trigger system, the Sum TriggerMJ. Recently, the VERITAS collaboration reported pulsed 
gamma rays at energies exceeding 100 GeV 151]. However, the cutoff energy of the Crab pulsar spectrum 
determined by Fermi-LAT was still ~ 6 GeV, which seemed to be in conflict with the very-high-energy 
data 1 6]. 

Here we present the spectral study of the Crab pulsar, using the public Fermi-LAT data and four years 
of MAGIC data recorded by the single telescope and the stereoscopic systems. 



2. The MAGIC telescopes 

The two MAGIC telescopes flU, situated on the island of La Pahna (28.8° N, 17.8° W, 2220ma.s.U, 
use the Imaging Atmospheric Cherenkov Technique to detect gamma rays above a few tens of GcVq 
During its monoscopic data taking era, a low-energy trigger system called Sum Trigger was developed 
and implemented in MAGIC-I in order to search for gamma-ray pulsars. It reduced the monoscopic 
energy threshold down to 25 GeV. 

Since MAGIC started operating in stereoscopic mode in summer 2009, its background suppression 
was substantially improved and a sensitivitjH of 0.8% Crab nebula units above 250 GeV has been 
achieved l^j. Since for technical reasons, the Sum Trigger cannot participate in the stereo trigger, 
MAGIC has two independent observation modes for pulsars: Monoscopic, with lower threshold, and 
stereoscopic, with lower systematics and better sensitivity. 



3. Mono-mode observations 

MAGIC-I observed the Crab pulsar with the Sum Trigger in the winters of 2007/08 and 2008/09. After 
a careful data selection, the total effective observation time was 59 hours. 

The light curve of the Crab pulsar obtained with the mono-mode observations is shown in Fig. [T] 
(left). Following the usual convention adopted in the EGRET era of PIr [0.94 - 0.04] and P2e 
[0.32 - 0.43], the numbers of excess events in PIe and P2e are 6200 ± 1400 (4.3 cr) and 11300 ± 1500 
(7.4cr). Summing up PI and P2, the excess corresponds to a significance of 7.5cr. The background level 
was estimated using the phase interval [0.52 - 0.88]. 

Based on these excess events, the phase-resolved energy spectra of the Crab pulsar above 25 GeV 
were computed and are shown as yellow squares in Fig. [21 They are compatible with power laws [11]. 
The energy spectrum measured by Fermi-LAT is also shown in the same figure. For the Fermi-LAT 
points, 1 year of the Fermi-LAT data (August 2008 - August 2009) were used. The MAGIC-Mono 
measurements are not compatible with an extrapolated exponential cutoff spectrum as determined by 
Fermi-LAT. A detailed statistical analysis showed the inconsistency amounts to 6.7cr, 3. Oct, and 5.8cr for 
(PI + P2)e, PIe and P2e, respectively (see lll2ll and 111 ill for more details on the monoscopic analysis). 



4. Stereo-mode observations 

In stereoscopic observations, MAGIC has collected 73 hours of optimal quality Crab pulsar data between 
October 2009 and February 2011. The right panel of Fig. [T] shows the light curve of the Crab pulsar 
obtained by stereo-mode observations in the total energy range of 46 - 416 GeV, and two energy sub- 
bins of it. The significance of the pulsation was tested with the Z^q test, the H test and a simple 
X^-test None of these makes an a-priori assumption concerning the position and shape of the pulsed 
emission, and they yield significances of 8.6 cr, 6.4 o" and 7.7 cr, respectively. By fitting two Gaussians 
to the two peaks, the peak positions are estimated to be 0.005 ± 0.003 and 0.3996 + 0.0014, while the 
corresponding FWHMs are 0.025 ± 0.007 and 0.026 ± 0.004. Defining the signal phases as ±2cr of the 
fitted Gaussian around the peaks, we calculate spectra both for the conventional, unbiased P1/2e intervals 
(see above), and for the narrower peaks PIm [0.983 - 0.026] and P2M [0.377 - 0.422] that we find in our 
lightcurve. 

The phase-resolved energy spectra are also shown in Fig. |2l They connect with the mono-mode 
measurements within statistical and systematic errors, and are not incompatible with a power law. Further 
details of the results of the stereo-mode observations are presented in il4il . 

' The threshold in standard trigger mode, defined as the peak of the simulated energy distribution for a Crab-nebula-like 
spectrum after all cuts and at low zenith angles, is 75 - 80 GeV. 

^ Defined as the source strength needed to achieve A'ex/ -^A'bkg = 5 in 50 h effective on-time. 
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Figure 1. MAGIC folded light curves of the Crab pulsar. Left: Monoscopic hght curve (2007 - 2009). 
Shaded areas are the on-phase regions P1/2e. Right: Stereoscopic light curve (2009 - 2011) for the total 
range in estimated energy, and for two sub-bins separately. The shaded areas are the on-phase regions 
PIm and P2m (see text), the light shaded area is the off-region [0.52 - 0.87]. The dashed Une is the 
constant background level calculated from that off-region. 



5. Discussion and conclusion 

We found a pulsed VHE gamma-ray signal from the Crab pulsar in monoscopic and stereoscopic MAGIC 
data that allows us to present spectra with an unprecedentedly broad energy range and phase resolution. 
Our spectra range over more than one order of magnitude, and, along with Fermi-LAT data [6], comprise 
for the first time a gamma-ray spectrum of the Crab pulsar from 100 MeV to 400 GeV without any gap. 
On the high-energy end, they are in good agreement with recently published VERITAS spectra of P14-P2 
above 100 GeV [5], including also the positions and narrow widths of the two pulses. 

A possible theoretical explanation of the spectrum iflll [l4ll is that the spectral tails arise from inverse 
Compton scattering of secondary and tertiary electron pairs on magnetospheric IR-UV photons. The 
model predicts the onset of that component at a few tens of GeV, producing a power law spectrum that 
possibly extends up to TeV energies. The modeling from il4\ is shown as a blueish solid line in Fig. [2a). 
It is compatible with the measurements of mono- and stereoscopic MAGIC, and VERITAS. It shall be 
noted that the GeV flux of the model in il4|] is too high to match the Fermi-LAT data because it contains 
bridge emission. 

The two crucial spectral features to establish in order to benchmark this and other models are a 
possible upward-kink in the transition region between curvature and hard component, and the detection 
or exclusion of a terminal cutoff at a few hundred GeV. The MAGIC telescopes, which are being upgraded 
in 2011, can provide both of these benchmarks in the coming years, when more data will improve the 
statistical precision of the measurements. 



6. Acknowledgments 

We would like to thank the Instituto de Astrofisica de Canarias for the excellent working conditions 
at the Observatorio del Roque de los Muchachos in La Palma. The support of the German BMBF 
and MPG, the Italian INFN, the Swiss National Fund SNF, and the Spanish MICINN is gratefully 
acknowledged. This work was also supported by the Marie Curie program, by the CPAN CSD2007- 
00042 and MultiDark CSD2009-00064 projects of the Spanish Consolider-Ingenio 2010 programme, by 
grant DO02-353 of the Bulgarian NSF, by grant 127740 of the Academy of Finland, by the YIP of the 



(a) Crab Pulsar, P1+P2 



ti 



(P1+P2 
(P1 + P2, 
(P1+P2, 
(P1+P2: 

(P1+P2: 
(P1+P2: 



, MAGIC stereo I 
MAGIC Stereo ( 
MAGIC Mono (A 

, Feimi-LAT (Abe 



etal. 2012] 
!tal. 20121 
: al. 20111 
2010) 
!. 20111 



ni-LAT (Ail 
Iv VERITAS (Aliu 
MAGIC Stereo (, 
MAGIC Stereo, ICRC2011 
Fermi-LAT (Abdo et al. 2010| 
ed, 0G+pair3 (Alek5i£ et al. 2012; 



I. 20121 



10^ 10^ 
Energy [ GeV ] 



(b) Crab Pulsar, PI 



ttl 



(c) Crab Pulsar, P2 



Figure 2. Compilation of spectral measurements of MAGIC and Fermi-LAT for the two emission peaks 
PI and P2 separately, and both peaks together. The VERITAS spectrum is only available (and shown) for 
Plv+P2v, which are yet narrower than P1/2m For comparison, the Crab nebula measurements of MAGIC 
and Fermi-LAT are also shown. Points of similar color refer to the same phase intervals (marked by the 
indices MAGIC, EGRET, VERITAS, see text). The MAGIC flux points are bias-corrected using an 
unfolding technique jlSll . The blueish soUd line is the modelling discussed in Section |5] 
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